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Exact calculations of the transmittance of surface corrugated optical waveguides are presented. 
The elastic scattering of diffuse light or other electromagnetic waves from a rough surface induces a 
diffusive transport along the waveguide axis. As the length of the corrugated part of the waveguide 
increases, a transition from the diffusive to the localized regime is observed. This involves an analogy 
with electron conduction in nanowires, and hence, a concept analogous to that of "resistance" can be 
introduced. We show an oscillatory behavior of both the elastic mean free path and the localization 
length versus the wavelength. 
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The influence of inhomogeneities and defects in the 
transmission properties of waveguides for both quantum 
and classical waves is a subject of increasing researcHJTj. 
The consequences of volume defects on electron wave 
propagation through two .dimensional (2D) wires have 
been recently addressecuJcl On the other hand, some 
aspects of the effects of slight surface roughness and vol- 
ume inhomogeneities both in optical fibers, and waveg- 
uides have been under study for a long timePH 11 !, mainly 
in connection with their attenuation; however, the study 
of these effects has been to date concentrated on the wave 
scattering into the radiation field; therefore, no analysis 
has been made so fap-pf the coupling between forward 
and backward modesEHl, neither on the influence of sur- 
face roughness induced localization effectstU into both 
this coupling and the radiation losses. 

Although some analogies between electron waves and 
light andx>ther classical waves have been put forward in 
aperturesH, in this letter we address localization effects 
of light and other electromagnetic waves propagating 
through a 2D corrugated waveguide (see inset of Fig. 1). 
There are two main differences between electron trans- 
port and the common situation in optical guides. First, 
the electron waves are confined within the conductor and, 
for small voltages, they do not couple with the radiation 
field. Second, the electron modes are mutually incoher- 
ent. So to resemble this picture we consider a waveguide 
with perfectly conducting walls, and the incident modes 
to be mutually incoherent. This can be achieved for ex- 
ample by previously passing the incident light through 
a moving diffuser so that correlation between phases of 
different modes is lost. 

The corrugated part of the waveguide, of total length 
L, is composed of n slices of length I. The width of 
each slice has random values uniformly distributed be- 



tween Wo — S and Wq + S about a mean value Wo- We 
shall take W /S = 7 and 1/6 = 3/2 (i.e. 25/W = 0.286 
and I /Wo = 0.214). The main transport properties do 
not depend on the particular choice of these parameters, 
however, s-polarization with the electric vector paral- 
lel to the grooves (TE modes) is assumed. Transmission 
and reflection coefficients are exactly calculated on solv- 
ing the 2D wave equation by mode matching at each slice, 
together with a generalized scattering-matrix technique, 
whose details can be found in Ref. 14. 

We define the normalized transmittance of the waveg- 
uide as G = Yljj Tij , where Ty is the ratio between the 
total flux transmitted from the incoming mode i into the 
outgoing mode j, and the total flux of this incoming 
mode, $™, i.e. = $° ut /$l n , and the sum runs over 
the total number of propagating modes. We also define 
the optical analogue of the "resistance" as the inverse of 
the transmittanceEl R — l/G. Ensemble averages, de- 
noted by (.), are performed over 100 realizations of the 
corrugated waveguide. 

For a perfect waveguide of constant width Wo, the nor- 
malized transmittance (G) is simply given by the number 
of propagating modes N, i.e. (G) = N = E{2Wo/A}. A 
plot of (G) versus Wo /A (see Fig.l) shows a staircase be- 
havior which constitutes the analogous for classical waves 
of the conductance quantization of electronic systemsta. 
Fig. 1 also contains (G) versus Wo/X for different length 
values L of the corrugated portion of the waveguide. For 
moderate lengths L, the transmittance shows a dip just 
at the onset of a new propagating mode. Similar |r dips 
appear in the case of scattering by volume defectsEO. 

As the length L of the corrugated region increases, (G) 
presents an oscillatory behavior that, as we show below, 
reflects successive transitions from diffusive to localized 
regimes. This is done on analysing the dependence of 
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(R) and (ln(G)) on L. Note that, while in clean sys- 
tems the transport is ballistic and consists essentially of 
unscattered waves, (the mean free path £ is then much 
larger than both L and Wo), in a randomly corrugated 
waveguide the two regimes, aforementioned above, are 
clearly visible as shown in Fig. 2, which contains (R) and 
(ln(G)) versus L, (cf. Fig. 2(a) and 2(b), respectively), 
for Wo/ A = 2.6. As L increases from zero, (R) follows 
first an almost perfect linearJpehavior with L which cor- 
responds to the semi-ballisticB regime. In this case, there 
is so much scattering in the longitudinal direction that 
the transport is diffusive while in the transversal direc- 
tion almost no scattering occurs (Wo < £ < L). Then, 
in analogy with Ohm's law for electron wires, the waveg- 
uide "resistance" can be described by an "ohm-like" tw#i 
proportional to L/Wq plus a "contact resistance" R c 



(i) 



where thei-i'resistivity" , p, and the mean free path £ are 
related byH £/W Q = l/(Np). 

This linear behavior breaks down as L increases fur- 
ther. Then, a typical linear decrease of (ln(G)) with L 
appears, as shown in Fig. 2(b). The system has now en- 
tered in the localized regime at which .Wp < £ <C L. The 
localization length £ is defined byOEj 



d(HG)) \ 
dL ) 



(2) 



The contact resistance R c and the effective mean free 
path £/Wq can be easily obtained both from Eq. [I] and 
by a least-square fitting of the linear part in Fig. 2(a). 
In Fig. 3(a) we show these values of R c and £/Wq versus 
Wo/A. The contact resistance (see the inset of Fig. 3) 
is slightly higher than the expected value 1/N, due to 
the transition region between ballistic and semi-ballistic 
regimes. The mean free path oscillates as Wo / A increases, 
having its minima close to half-integers of Wo/A, which 
correspond to the appearence of a new propagating mode 
in the waveguide, and thus establishes the above quoted 
connection with the transmittance dips at the onset of 
new modes shown in Fig. 1. 

The dependence of the localization length £ on the 
wavelength (Fig. 3(b)) can be also obtained from least- 
square fitting of the linear part of the plot of (ln(G)) 
versus L in Fig. 2(b). It should be noted that, within 
the numerical accuracy of our calculation, £ « N£ as 
shown in Fig. 3(b). This again has an analogy in agree- 
ment with the expected beha vior for electron transport 
through disordered mediallal23. 

The wavelength dependence of both the localization 
length and the mean free path shows that, for a fixed 
length, it is possible to pass on from the localized regime 
into the diffusive one by changing the wavelength of the 
incoming wave. 

In conclusion, we have shown that diffuse light and 
other electromagnetic wave propagation through surface 



corrugated perfectly reflecting waveguides, constitutes 
the optical analogue of electron transport in disordered 
nanowires. As such, the concept of "averaged resistance" 
and the optical analogous of "Ohmic" behavior in the 
semi-ballistic (diffusive) regime can be introduced for 
these systems. We have also shown that both the effec- 
tive mean free path and the localization length oscillate 
with the wavelength. This oscillatory transition from dif- 
fusive to localized behavior should be observed in optical 
waveguides. The influence of surface induced localizaton 
effects on radiation losses in actual waveguides remains 
an open question. We hope that our results will estim- 
ulate both experimental and theoretical research in this 
direction. 
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FIG. 1. Averaged transmitance (G) as a function of Wo/\ 
for fixed values of 5/W and l/Wo, (Wo = 7 5 and I = (3/2) S). 
Note that the case L/Wo = corresponds to a flat waveguide. 
The inset shows a schematic view of the system. 

FIG. 2. (a) (R) versus the length L of the corrugated 
region. The straight line represents the best fitting of the 
linear behavior associated to the diffusive transport regime 
(see Eqjl|). The results are obtained for Wo/\ = 2.6, i.e they 
correspond to 5 propagating modes in the waveguide. The 
corresponding effective mean free path is l/Wo ~ 6.3. (b) 
(ln(G)} versus L for the same case as in (a). The best fitting of 
the linear part of the curve is also shown. The corresponding 
localization length is £/Wo 34.4. 

FIG. 3. (a) Effective mean free path l/Wo versus Wo/\. 
The inset shows the behavior of the contact resistance R c to- 
gether with the expected value for a perfect waveguide (dot- 
ted line), (b) Localization length £/Wo (filled circles) and Nl 
(open circles) versus Wo/X. £ ~ Nl, within the numerical 
accuracy. 



3 




0.5 1.0 1.5 2.0 



Fig. 1 
A. Garcfa-Martm, et al. 




Fig. 3 
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